Journal  of  Power  Sources  185  (2008)  1141-1146 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Journal  of  Power  Sources 

journal  homepage:  www.elsevier.com/locate/jpowsour 


pri 

SbbaautS 


Identification  problems  and  analysis  of  the  limit  current  in  fuel  cells 

Paolo  Costa,  Barbara  Bosio* 

Department  of  Civil,  Environmental  and  Architectural  Engineering,  University  of  Genoa,  Via  Opera  Pia  15,  16145  Genoa,  Italy 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  30  March  2008 
Received  in  revised  form  9  July  2008 
Accepted  13  July  2008 
Available  online  25  July  2008 


Keywords: 

Mass  transfer 
Limit  current 
Utilisation  factor 
Fuel  cell  performance 


This  paper  presents  our  analysis  of  the  mass  transfer  phenomena  in  fuel  cells.  Our  attention  was  specif¬ 
ically  focused  on  the  parameter  identification  problems  (briefly  identification  problems)  related  to  the 
interpretation  of  the  experimental  data  on  the  limit  current. 

Two  factors  work  together  to  make  the  interpretation  of  the  limit  current  data  difficult:  the  presence 
of  two  interacting  electrodes  (electrode  coupling  phenomena)  and  the  flow  rate  changes  along  each  elec¬ 
trode.  The  first  effect  can  be  well  defined  and  understood  by  studying  reagent-poor  inlet  flow  rates  to  the 
electrodes;  the  second  effect  can  be  analysed  by  superimposition.  In  this  analysis,  the  asymptotic  solu¬ 
tions  for  complete  decoupling  conditions,  that  are  complete  anodic  (or  cathodic)  mass  transfer  control, 
have  been  shown  to  be  very  useful  reference  tools. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Today  fuel  cells  are  in  a  critical  development  stage,  which  can 
be  facilitated  with  a  work  aimed  at  investigating  basic  phenom¬ 
ena  more  deeply  as  well  as  improving  technologies  and  processes 
optimising  performance,  lifetime  and  costs. 

In  this  scenario,  a  careful  study  of  the  mass  transfer  resistances 
can  be  interesting  because  operation  under  limiting-diffusion  con¬ 
ditions  can  strongly  penalise  fuel  cell  functioning. 

Usually  studies  focused  on  the  improvement  of  fuel  cell  perfor¬ 
mance  take  account  of  the  effects  of  the  mass  transfer  resistances  in 
terms  of  the  concentration  polarisation,  the  fuel  utilisation  factor 
or  the  limit  current  [1-3],  nevertheless  the  subject  is  not  usually 
approached  from  the  point  of  view  of  experimental  data  interpre¬ 
tation.  In  this  paper  such  an  approach  will  be  discussed. 

Limit  current  experimental  data,  although  relatively  rare, 
are  the  most  specific  and  direct  way  to  characterise  the  mass 
transfer  resistances  of  a  cell,  because,  by  definition,  the  limit 
currents  are  completely  controlled  by  these  resistances,  while 
they  are  completely  independent  of  the  intrinsic  electrochemical 
kinetics. 

Mass  transfer  resistances,  however,  cannot  be  automatically  and 
directly  determined  by  limit  current  data  without  using  some  inter¬ 
pretative  procedure.  In  fact,  two  factors  work  together  to  make  the 
problem  of  mass  transfer  identification  difficult:  the  electrode  cou¬ 
pling  phenomena  and  the  flow  rate  changes  along  each  electrode. 
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The  first  effect  corresponds  to  a  mixed  mass  transfer  control;  it 
acts  when  only  a  portion  of  the  cell  is  controlled  by  anodic  mass 
transfer  phenomena  (e.g.  H2  is  the  limiting  reactant)  while  the 
remaining  portion  is  controlled  by  cathodic  ones  (e.g.  oxidant  is 
the  limiting  reactant)  [4].  This  effect  can  be  isolated  from  the  sec¬ 
ond  by  studying  reagent-poor  inlet  flow  rates  to  the  electrodes 
and  determining  the  fraction  of  the  cell  surface  which  under¬ 
goes  anodic  (or  cathodic)  mass  transfer  control  under  such  simple 
conditions. 

The  second  effect  depends  on  the  non-linear  relationship 
between  the  utilisation  factor  and  the  reagent  concentration  in 
reagent-rich  systems:  as  already  shown  in  our  previous  works  [5,6] 
a  superficial  analysis  of  the  concentration  field  effects  can  lead  to 
significant  errors  in  the  estimation  of  the  mass  transfer  coefficients. 

In  the  following  sections  these  two  effects  and  their  super¬ 
imposition  will  be  discussed  in  order  to  provide  correct  data 
interpretation. 

The  work  refers  in  particular  to  molten  carbonate  fuel  cell 
(MCFC)  technology,  but  its  matter  is  valid  for  each  kind  of  fuel  cells. 

2.  The  reference  case 

Reference  has  been  made  to  a  cross-flow  fuel  cell  (in  partic¬ 
ular  an  MCFC,  which  is  a  molten  carbonate  fuel  cell)  of  a  typical 
industrial  laboratory  scale  (with  a  rectangular  Lx  x  Ly  shape,  about 
0.1  m  x  0.1  m)  [6].  Such  laboratory  devices  can  be  considered  sub¬ 
stantially  isothermal  and  isobaric,  while  their  flow  field  can  be 
schematised  as  a  cross-coupling  of  two  plug-flow  systems,  namely 
an  anodic  one  flowing  from  left  to  right  and  a  cathodic  one  flowing 
upwards. 


1142 


P.  Costa,  B.  Bosio  /  Journal  of  Power  Sources  185  (2008)  1141-1146 


Nomenclature 

C  reactant  concentration  (mol  m-3 ) 

/A  see  Eq.  (19) 

fc  see  Eq.  (20) 

F  Faraday’s  constant  (A  s  mol-1 ) 

kc  mass  transport  coefficient  (m  s-1 ) 

/  cell  electrical  current  (A) 

J  cell  electrical  current  density  (A  m-2 ) 

Lx  length  of  cell  side  parallel  to  anodic  flow  (m) 

Ly  length  of  cell  side  parallel  to  cathodic  flow  (m) 

n  number  of  electrons 

q  volumetric  flow  rate  (m3  s-1 ) 

r  reaction  rate  per  unit  surface  (mol  m-2  s) 

u  utilisation  factor 

x  dimensionless  coordinate  of  the  anodic  flow,  Eq.  (4) 

x!  see  Eq.  (24) 

X  coordinate  of  the  anodic  flow  (m) 

y  dimensionless  coordinate  of  the  cathodic  flow,  Eq. 

(4) 

y'  see  Eq.  (24) 

Y  coordinate  of  the  cathodic  flow  (m) 

z  molar  fraction 

Greek  letters 

f  see  Eq.  (4)  (mol  m-2  s) 

vA  see  Eq.  (18) 

vAi  stoichiometric  coefficient  of  component  i  at  the 

anode 

v ala  stoichiometric  coefficient  of  limiting  reactant  at  the 

anode 

Vq  see  Eq.  (20) 

vci  stoichiometric  coefficient  of  component  i  at  the 
cathode 

v clc  stoichiometric  coefficient  of  limiting  reactant  at  the 

cathode  e 

cpc  fraction  of  the  cell  surface  under  cathodic  control 

@c  fraction  of  the  cell  surface  under  cathodic  control 

O'  see  Eq.  (A6) 

Subscripts  and  superscripts 
A  anode 

b  boundary 

C  cathode 

LA  limiting  reactant  at  the  anode 

LC  limiting  reactant  at  the  cathode 

o  inlet 

t  total 

intersection  point  of  asymptotic  solutions  (12)  and 
(13) 


Moreover,  for  a  first,  simplified  approach,  reagent-poor  solu¬ 
tions  and  invariant  volumetric  flow  rates  can  be  assumed  for  the 
anodic  (uniform  qA)  and  the  cathodic  (uniform  qc)  sides,  respec¬ 
tively.  Finally,  explicit  reference  has  been  made  to  the  limit  current 
condition,  which  has  been  assumed  over  the  entire  cell  plane,  so 
that  the  local  kinetics  are  everywhere  controlled  by  transport  phe¬ 
nomena  and,  more  precisely,  by  which  of  the  two  reactant  diffusive 
transport  mechanisms  is  slower,  the  anodic  or  the  cathodic: 

r  =  -4  =  min(fccLACLA,  kcLCCLC)  (1 ) 


For  instance,  in  an  MCFC,  where  the  reactions 
H2  +  C032-  H20  +  C02 +  2e_  (at  the  anode) 

C02  +  \02  +  2e-  C032-  (at  the  cathode) 

occur  (neglecting  the  possibility  of  direct  oxidation  of  CO  or  hydro¬ 
carbons),  the  anodic  reactant  is  hydrogen  (LA  =  H2 )  and  the  cathodic 
limiting  reactant  is  normally  carbon  dioxide  (LC  =  C02),  while  the 
second  cathodic  reactant,  oxygen,  is  present  in  excess  and  does  not 
show  any  limiting  effect. 

In  such  simple  circumstances  the  local  mass  balances  of  the 
anodic  reagent  LA  and  the  cathodic  reagent  LC  can  be  written, 
respectively  as 

dC 

=  _m^n^<cLACLA’  kdcQ x)h>  0  <X  <  Lx, 

X  =  0;  CLA  =  CLAo  (2) 

and 

Qc~^  =  -min  (/<cLACLA,  kcLCCic)Lx,  0  <Y  <  Ly, 

Y  =  0;  CLC  =  CLCo  (3) 

or  in  the  compact  forms 

f A  =  kcLaCLA  ’  X  =  ^cLA  ~~~  •>  xt  =  W<cLA  ~~~  \ 

Qa  Qa 

fc  =  /<clcClc  y  =  Yk&c  fri  yt  =  Lykc lc  ~r~  (4) 


d?A 

dx 


Ly 

-min(fA,fc),  0  <  x  <  xt  =  kc  la^x ~~ >  x  =  0;  ^A  =  £ao 

Qa 

(5) 


=  -min(fA,  fc),  0<y<yt  =  kCLc£yA 
oy  qc 

y  =  0;  fc  =  fco  =  ^ao  (6) 

introducing  the  parameter  ^  =  fco/fAo- 

3.  Simple  and  mixed  control  conditions 

According  to  Eqs.  (5)  and  (6),  the  plane  of  the  cell  can  be 
divided  into  two  zones  (Fig.  1):  one,  where  ^ak^o  which  is  con¬ 
trolled  by  the  anodic  mass  transfer  phenomena,  and  one,  where 
fc  <  f a.  which  is  controlled  by  the  cathodic  transfer  phenomena. 
The  boundary  between  these  two  zones  is  obviously  defined  by  the 
condition  f  A  =  f c  and,  then,  described  by  the  straight  line  equation 
(see  Appendix  A): 


y  =  n  1+X)-1,  *  =  i<clcClco  (7) 

SAo  '^cLAClAo 

or,  in  terms  of  the  intercept  of  the  boundary  on  the  left  margin  of 
the  cell: 

y  =  yx=o  +  y*=o  =  v- 1  (8) 

The  comparisons  of  theyx=0  value  given  by  the  second  Eq.  in  (8) 
and  the  yt  value  given  by  the  second  Eq.  in  (6)  and  of  the  y  value 
given  by  the  first  Eq.  in  (7)  and  the  xt  value  given  by  the  second 
Eq.  in  (5)  show  that  the  cell  can  undergo  three  different  kinds  of 
current  limitations: 
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Fuel 


vp  parameter  [-] 


Fig.  1.  Cross-flow  fuel  cell  under  limit  current  conditions,  (a)  Changing  of  the  bound¬ 
ary  position  for  different  values  of  parameter  tfC  (b)  Cell-surface  fraction  working 
under  cathodic  control  conditions  versus  parameter  ^  (underlined  as  an  example 
the  boundary  for  V  =  0.75).  For  a  cell  with  yt  =  1,  xt  =  2,  the  complete  anodic  control 
requires  ^  >  2,  while  the  complete  cathodic  control  requires  ^  <  1/3.  In  intermediate 
conditions  (1/3  <  ^  <  2),  the  cell  works  under  mixed  control. 


•  Simple  anodic  mass  transfer  control: 

for  &  >  1  +  yt  (9) 

The  boundary  position  is  characterised  by  y  >yt  for  0  <  x  <  xt :  the 
boundary  lies  to  the  left  of  the  upper-left  edge  of  the  cell  and  does 
not  cross  it,  so  that  £A  <  £c  and  the  anodic  mass  transfer  controls 
each  point  of  the  cell.  In  such  case  the  boundary  can  be  called 
“external”. 

•  Simple  cathodic  mass  transfer  control: 

for  tf'c  (1 +xt)-1  (10) 

The  “external”  boundary  position  is  now  characterised  by  y  <  0 
for  0  <  x  <  xt:  the  boundary  lies  to  the  right  of  the  lower- right  edge 
of  the  cell  and  again  does  not  cross  it.  In  this  case  £c  <  ?a  and  the 
cathodic  mass  transfer  controls  each  point  of  the  cell. 

•  Mixed  mass  transfer  control: 

for  (1  +xt)“1  <  V<  1  +yt  (11) 

In  this  case  the  boundary  crosses  really  the  cell  and  divides 
it  into  two  parts,  the  upper  under  anodic  control  and  the  lower 
under  cathodic  control.  Under  such  conditions,  the  fraction  of  the 
cell  surface  working  under  cathodic  control  increases  from  zero 
to  unity  with  the  decrease  in  parameter  &  in  the  range  of  Eq.  (11 ). 
An  example  of  this  trend  is  given  in  Fig.  1 ,  which  refers  to  a  cell 
“shape”  xt  =  1,  yt  =  2  and  a  mixed  control  range  1/3  <  ^  <  2. 


Fig.  2.  Cross-flow  fuel  cell  under  limit  current  conditions:  an  example  of  the  trend 
of  the  ratio  £b/£Ao  along  the  boundary  between  the  anodically  and  cathodically  con¬ 
trolled  zones,  foryt  =  1,  xt  =  2  and  different  values  of  &  ( &  =  1.5  for  0  <  x  <  0.33 ;  ^  =  1.25 
for0<x<0.6;  «^  =  lfor0<x<l;  ^  =  0.75 for  0.33  <x<  1.66;  ^  =  0.5 for  1  <x<2;  «^  =  0.4 
for  1.5  <  x  <  2).  When  ^  >  2  (complete  anodic  control)  and  ^  <  1  /3  (complete  cathodic 
control)  the  boundary  does  not  cross  the  cell. 


4.  The  limit  performance  of  the  cell 


The  performance  of  a  cell  working  under  a  simple  mass  transfer 
control  can  be  straightforwardly  predicted:  for  instance,  the  anodic 
utilisation  factor  for  a  cell  under  simple  anodic  control  is  [4] 

Ua  =  «la  =  1  -exp(-xt),  ^  >  1  +yt  (12) 


Obviously  a  similar  expression  holds  for  the  cathodic  utilisation 
under  simple  cathodic  control.  In  such  circumstances  the  anodic 
utilisation  factor  is 

Ua  =  «lc^,  Ulc  =  1 -exp(-yt),  ^<(l+xt)_1  (13) 

yt 


Under  mixed  control  conditions,  the  anodic  composition  lies  within 
the  limits  of  Eqs.  (12)  and  (13)  (see  Appendix  A).  An  example  of  the 
results  obtained  is  given  in  Fig.  3,  where  the  anodic  utilisation  is 
reported  as  a  function  of  the  parameter  &  and  the  limit  trends  for 
complete  anodic  control,  Eq.  (12),  and  complete  cathodic  control, 
Eq.  (13),  are  shown  as  well. 

In  Fig.  3  complete  anodic  or  cathodic  control  corresponds  to 
^>2  and  ^<1/3,  respectively,  but  it  is  evident  that  the  asymp¬ 
totic  solutions  (12)  and  (13)  can  be  considered  sufficiently  accurate 
over  wider  ranges  of  values,  namely  XI'>\  and  &  <  0.6.  Moreover, 
for  many  purposes  the  validity  of  the  asymptotic  behaviour  can  be 
extended  to  the  intersection  point. 


q,*  _  yt[l  -  exp(-xt)] 
xt[l  -  exp(-yt)] 


(14) 


/  yt  =  1,  Xj  =  2,  V* 

=  0.68 

/  V  <  0.33 

complete  cathodic  control 

/  *F>2 

complete  anodic  control 

/  0.33  <*F  <2 

mixed  control 

/  0.33  <^<0.68 

prevailing  cathodic  control 

/  0.68  <*P<2 

prevailing  anodic  control 

0  0,5  1  1,5  2  2,5 

lP  parameter  [-] 


Using  Eq.  (7)  the  value  of  £b(x)  acting  on  the  boundary  can  also 
be  obtained  (see  Appendix  A).  An  example  is  shown  in  Fig.  2. 


Fig.  3.  Cross-flow  fuel  cell  under  limit  current  conditions:  an  example  of  the  trend 
of  the  anodic  utilisation  factor  for  yt  =  1,  xt  =  2  and  for  different  values  of  flC 


1144 


P.  Costa,  B.  Bosio  /  Journal  of  Power  Sources  185  (2008)  1141-1146 


with  a  maximum  error,  just  at  the  intersection  point,  of  the  order 
of  10%.  In  such  a  way,  the  mixed  control  zone  can  be  divided  in  two 
parts: 

<  &  <  1  +yt,  uA  «  u la  (prevailing  anodic  control); 

(l+xt)-1  <  W  <  &*,  uA  %  (prevailing  cathodic  control) 

y  t 

(15) 


Another  way  of  arranging  the  above  results  is  to  relate  them  to 
a  cell  working  under  constant  inlet  conditions  in  terms  of  reactant 
concentration  for  one  electrode  (e.g.  the  cathode,  constant  CLCo) 
and  variable  inlet  conditions  for  the  other  (anode,  variable  CLAo). 

As  the  cell  current  is  proportional  to  the  anodic  utilisation  and 
the  anodic  reagent  inlet  rate 


I  =  nFqACLAouA  (16) 

the  limit  current  of  the  cell  can  be  obtained  by  Eq.  (16),  in  which  the 
limit  utilisation  is  expressed  by  equations  such  as  (A16),  (A17)  or  by 
the  approximations  (15),  (12),  (13).  In  particular,  for  constant  qA,  qc , 
CLco  and  variable  CLAo,  the  limit  current  will  be  proportional  to  uA 
and  the  inlet  ratio  Qa^laoIQc^ico  =  const./*/',  so  that  the  asymptotic 
trends  of  the  cell  will  be 

i !'>'!'*,  /  %  /A  =  I  </'<</'*,  I*IC  =  const.  (17) 

as  shown  in  Fig.  4. 

The  relevant  point  is  that  experimental  findings  are  char¬ 
acterised  by  this  kind  of  trend  in  terms  of  limit  current  and 
anodic  composition,  see  Fig.  5.  As  the  anodic  inlet  concentration  is 
increased,  the  limit  current  increases  too,  initially:  under  these  con¬ 
ditions  there  is  complete,  or  at  least  prevailing,  anodic  control,  in 
other  words  the  two  electrodes  are  completely  decoupled  and  only 
the  anode  is  significant.  After  a  narrow  mixed  control  range  (signif¬ 
icantly  coupled  electrodes),  the  limit  current  becomes  constant,  as 
it  is  controlled  by  the  constant  cathode:  now  the  cell  works  under 
complete,  or  at  least  prevailing,  cathodic  control  and  the  decou¬ 
pling  of  the  electrodes  is  again  good,  with  only  the  cathode  being 
significant  [4]. 

The  experimentally  confirmed  applicability  of  the  discussed 
analytical  approach  to  small-size  cells  puts  in  evidence  its  useful¬ 
ness  in  terms  of  analysis  and  diagnostic  aims. 

Considering  commercial  cells  of  full  area  (usually  about  1  m2), 
the  non-uniformity  of  the  chemico-physical  properties  on  the  cell 
plane  can  be  significant.  In  particular,  the  temperature,  flow  rate 
and  transport  coefficient  distributions  can  affect  the  cell  behaviour, 


Fig.  4.  Cross-flow  fuel  cell  under  limit  current  conditions:  the  cell  current  as  a 
function  of  the  anodic  inlet  under  constant  cathodic  conditions.  The  cell  current 
is  reduced  in  respect  to  the  complete  conversion  of  the  cathodic  flow. 


Fig.  5.  Cross-flow  fuel  cell  under  limit  current  conditions:  the  experimental  cell 
current  as  a  function  of  the  anodic  inlet  under  constant  cathodic  conditions. 


so  that  the  predictable  boundary  position  could  need  some  cor¬ 
rections.  On  the  other  hand,  it  is  interesting  to  observe  that  the 
boundary  linearity  is  well  confirmed  for  commercial  MCFCs  work¬ 
ing  under  standard  operating  conditions.  This  fact  has  been  verified 
by  means  of  simulation  with  a  detailed  model  [7]. 


5.  Flow  rate  effects 


As  electrode  reactions  are  not  equimolar,  the  molar  flow  rates 
can  undergo  significant  changes  in  flowing  along  the  electrode 
paths.  Under  constant  temperature  and  pressure  the  same  changes 
affect  the  volumetric  flow  rates  too.  For  instance,  the  anodic  flow 
rate  in  an  MCFC  increases  along  the  reaction  path  and  so  depends 
on  the  molar  fraction  of  the  hydrogen  present  [5].  In  terms  of  vol¬ 
umetric  flow  rate  and  concentration  the  dependence  is 


Qa  =  Qao 


C  +  fACLAo 


^  ^ALA 


C  +  vACLA  ’ 
or,  by  using  the  definition  of  tfA  in  Eq.  (4), 


Qa  _  1  +/a£ao  r  _  t>A 
Qao  1  +/a£a  ’  A  kcLAC 


(18) 


(19) 


It  is  evident  that,  as  vA  =  l  in  MCFCs,  for  a  dilute  reactant 
(CLa  ->  0)  the  volumetric  flow  rate  can  be  considered  as  a  constant, 
while  it  tends  to  a  100%  variation  if  pure  hydrogen  is  fed. 

Similarly,  the  volumetric  flow  rate  at  the  cathode,  where 
vc  =  -1.5  in  MCFCs,  will  decrease  according  to 

Qc  =  1  +/c£co  f  =  vc  =  vg/vclc  ,  . 

Qc o  1  +/c£c  ’  C  KicC  kcL cC 

so  that  assuming  now 

x  =  XkcLA-^,  y  =  Ykclc^~  (21) 

Qao  Qc O 

Eqs.  (5)  and  (6)  become 


(^a/9x)(1  +/a£ao) 
(1  +/a£a)2 
Oft/foXi  +Kco) 

(i  +/cCc)2 


-min(£A,  Cc), 


-min(£A,  £c)> 


x  6,  —  £ao 

y  =  o>  tc  =  Cco 


(22) 

(23) 
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Fig.  6.  Cross-flow  fuel  cell  under  limit  current  conditions:  an  example  of  the  effect 
of  changing  the  flow  rate  on  the  asymptotic  solutions  for  the  cell.  The  approximate 
asymptotic  solutions,  calculated  assuming  uniform  flow  rates  along  the  electrodes 
(dotted  lines),  are  compared  to  the  correct  ones  (continuous  lines). 


The  system  (22),  (23)  can  be  reduced  to  the  form  (5),  (6)  when 
the  differential  substitutions 


dx'  = 


(1  +/a?a)2  dx 
1  +/a£ao 


(1  +/c£c) 
(1  +/c£  Co) 


(24) 


are  introduced. 

In  such  a  way,  most  of  the  results  of  the  previous  section  can  be 
rewritten  and  extended  to  those  cases,  frequent  in  practice,  where 
the  hypothesis  of  uniform  flow  rates  has  to  be  rejected  in  a  fairly 
accurate  description.  In  particular,  the  anodic  utilisation  uA  or  the 
cell  current  I  can  be  calculated  for  any  feed  composition  (zLAo,  zLCo) 
when  the  transport  coefficients  (/<cLA,  /<clb)  are  known  or,  inversely, 
in  an  identification  problem,  a  best  fit  for  the  transport  coefficients 
can  be  performed  from  measured  I  data.  When  mixed  mass-transfer 
conditions  are  controlling,  however,  the  calculations  are  not  imme¬ 
diate  and  the  explicit  solutions  cannot  be  expressed  in  analytical 
form. 

In  particular,  the  asymptotic  solutions  (12)  and  (13)  are  substi¬ 
tuted  by  implicit  equations  such  as  the  following: 


•  Anodic  asymptotic  solution: 


6.  Conclusion 

In  this  paper  a  discussion  of  the  parameter  identification 
problems  affecting  a  mass  transport  phenomena  study  has  been 
presented  in  order  to  provide  a  guide  to  a  correct  analysis  of  the 
limit  current  in  fuel  cells.  A  particular  reference  to  molten  carbonate 
fuel  cells  has  been  done. 

Analysing  the  behaviour  of  a  laboratory  fuel  cell  under  limiting- 
diffusion  conditions  and  uniform  temperature,  a  boundary  on  the 
cell  plane  has  been  identified  distinguishing  the  area  under  the 
anodic  control  and  the  area  under  the  cathodic  control  and  so 
underlining  the  possibility  of  a  mixed  control.  This  boundary  has 
been  demonstrated  linear  when  the  flow-rate  variation  along  the 
cell  is  negligible. 

In  addition,  a  linear  trend  of  the  cell  limit  current  versus  reagent 
concentration  has  been  shown  when  complete  decoupling  condi¬ 
tions,  that  are  complete  anodic  (or  cathodic)  mass  transfer  control, 
are  considered.  Such  asymptotic  solutions  have  been  demonstrated 
a  useful  reference  tool,  which  can  involve  a  maximum  error  of  the 
order  of  10%  in  case  of  mixed  control. 

By  this  way  the  interpretation  of  the  limit-current  data,  which 
are  the  most  specific  and  direct  data  to  characterise  the  mass 
transfer  resistances  of  a  cell,  has  been  shown  complicated  by  two 
non-linearity  effects:  mixed  control  and  flow  rate  effects. 

The  interpretation  of  such  non-linearity  effects  has  been  dis¬ 
cussed  and  shown  to  be  difficult  to  interpret,  as  the  two  effects 
tend  to  obscure  each  other. 

In  the  light  of  such  considerations,  detailed  local  studies  and  spe¬ 
cific  tests  involving  complete  electrode  phenomena  decoupling  and 
reagent-poor  feeds  are  advisable  for  evaluating  diffusion  contribu¬ 
tions.  Using  this  approach  it  is  possible  to  obtain  a  better  analysis 
of  the  fuel  cell’s  performance  and,  consequently,  more  effective 
performance-improvement  policies  can  be  set  up. 

Appendix  A 

A  A.  The  position  of  the  boundary 

The  boundary  is  defined  by  the  position  £A  =  £b  =  fb-  Moreover, 
starting  from  the  point  (x,  0)  and  following  the  cathodic  flow 
upward  (x  =  const.)  through  the  anodic  control  zone  to  the  boundary 
point  (x,  yb),  we  have 


ua  =  Ula,  jjp  =  QaoClaoUla 

with  xt  =  -[zlaoVaUla  +  (1  +  zLAovA)  ln(l  -  14*)] 


•  Cathodic  asymptotic  solution: 


^LcgcoCLCo 

QAo^LAo 


=  QCoClCoUic 


(25) 


ryb 

?b  =  £bo  -  /  - 

Jo 


?a  d y 


(Al) 


and,  by  differentiation  and  incorporating  Eq.  (5)  in  the  above  Eq. 
(Al), 


/ 


_  ^  ^  _  £A<tyb 

dxb  /„  dx  dxb 


l 


^Ady-^r  =  ^c°-^ 


with  yt  =  -[ZlcoI’cI'lc  +  (1  +zLc0vc)  ln(l  -  uLC)]  (26) 


fbdyb 

dxb 


(A2) 


In  such  a  case,  if  a  cell  is  working  under  complete  anodic  con¬ 
trol,  the  dependence  of  the  limit  current  on  the  composition  of 
the  anodic  feed  is  no  longer  a  simple  proportionality  (see  Fig.  4). 
On  the  contrary,  significant  non-linear  effects  are  evident  as  the 
hydrogen  content  increases:  an  example  is  shown  in  Fig.  6,  where 
the  approximate  asymptotic  solutions,  calculated  assuming  uni¬ 
form  flow  rates  along  the  electrodes  (dotted  lines),  are  compared 
with  the  correct  ones  (continuous  lines).  The  example  refers  to  a 
constant  inlet  composition  at  the  cathode  (ZLCo  =  0.2)  and  inlet  flow 
rate  ratio  (qco/^Ao  =  5),  the  values  being  chosen  with  reference  to  a 
typical  run  of  an  laboratory  MCFC  [6]. 


Similarly,  starting  from  the  point  (0,y)  and  following  the  anodic 
flow  from  left  to  right  (y  =  const.)  into  the  cathodic  control  zone  to 
the  boundary  point  (xb,  y),  we  have 


d?b  _  f  ,  <bdxb 

d^-?A°-fb-^T 


(A3) 


By  combining  Eqs.  (A2)  and  (A3),  it  is  possible  to  obtain  the 
differential  equation  of  the  boundary  line: 


dyb  _  £co 

dxb  £a° 


(A4) 
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The  boundary  condition  for  Eq.  (A4)  can  be  obtained  from  Eq. 
(Al)  written  for  x  =  0: 

x  =  o,  yb  =  y*=o,  £ao  —  £co —  ^Aoyx=o 

By  putting 

i p  — 

Kao 

Eqs.  (A4)  and  (A5)  can  be  rewritten 

Z-*'  x-° ■ 

so  that,  by  integration,  the  boundary  can  be  described  by  the 
straight  line  equation: 


(A7) 


yb  =  +  V  -  1 

A.2.  The  degree  of  cathodic  control 


(A8) 


Starting  from  Eq.  (8),  elementary  geometrical  considerations 
make  it  possible  to  evaluate  the  fraction  cpc  of  the  cell  surface  which 
is  under  cathodic  control.  For  instance,  under  the  condition  xt  >yt, 
we  have  (see  Fig.  1) 

(pc  =  0,  for  &  >  1  +yt 


(pc  = 
(pc  = 


(yt-^  +  ir 

2^ytxt 
yt-2^  +  2 


for  1  <  ^  <  (1  +yt) 


,  for  (1  +yt)(l  +xt)  1  <  V  <  1 
2 


2Vxt 


(Pc  =  ,  for  (1  +xt)  1  <  ^<  (1  +  yt)(l  +Xt)  1 


2  ^ytxt 
(pc  =  1,  for  ^  <  (1  +xt) 


-l 


(A9) 


A.3.  The  composition  on  the  boundary 

By  combining  the  differential  Eqs.  (A2)  and  (A4)  with  the  posi¬ 
tion  (A6),  the  changing  composition  along  the  boundary  can  be 
described  as 


d£b 

dxb 


=  Vito  -  ?b(l  -  V) 


(A10) 


The  starting  condition  for  Eq.  (A10)  can  be 


The  integration  of  the  Eq.  (A10)  with  the  condition  (All )  or  (A12) 
gives,  respectively: 


(A5) 

ft 

fto 

=  ^  +  exp[-0  +  ^)xbl)  1<;^<1+yt 

(A13) 

(A6) 

ft 

&  +  W2  exp[(l  -  'I'2)/'!']  exp[-(l  +  &)  xb] 

^Ao 

i + <^ 

(1  +xt)_1  <V<\ 

(A14) 

A.4.  The  performance  of  the  cell  in  terms  of  utilisation 

Under  simple  control  conditions  the  anodic  utilisation  factor  of 
the  cell  is  well  described  by  Eqs.  (12)  and  (13),  for  complete  anodic 
and  cathodic  controls,  respectively.  Under  mixed  control  condi¬ 
tions,  the  anodic  composition  at  a  point  on  the  right  side,  which 
is  delimited  by  the  boundary  and  where  the  anodic  mass  transfer 
is  controlling,  is 

ft(x)  =  fto  exp(-x),  xb  <  0;  ?A(x)  =  &  exp(xb  -  x),  xb  >  0 

(A15) 

so  that  the  mean  outlet  composition  is 


ftm  =  W1  -  UA)  =  /  ft(*t)  dy 


T 


(A16) 


Substituting  Eqs.  (A13)  or  (A14)  with  fb  and  making  the  integral 
in  Eq.  (A16),  the  following  results  can  be  obtained: 

1  -uA  =  exp(-x),  ^  >  1  +yt 

exp(-x){<^2  exp[(yt/^)-l  +  l/^]-exp(^-l-yt)) 

1_Ua  = - OTWt - ’ 

1  <  ^  <  1  +  y  t 

exp(-x){<^2  exp[(yt/</')  - 1  +  l/^]-exp[(^/l)-l  —  yt)]} 


1— uA=- 
(1  +yt)/(l  +xt)  <  V<  1  +yt 


(1  +<#0yt 


(A17) 


and  so  on,  from  elementary  symmetry  considerations,  as  ^ 
decreases  to  (1  +  xt)_1. 

An  example  of  the  results  obtained  is  given  in  Fig.  3. 
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